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There can exist topological obstructions to continuously deforming a gapped Hamiltonian for free 
fermions into a trivial form without closing the gap. These topological obstructions are closely 
related to obstructions to the existence of exponentially localized Wannier functions. We show that 
by taking two copies of a gapped, free fermionic system with complex conjugate Hamiltonians, it is 
always possible to overcome these obstructions. This allows us to write the ground state in matrix 
product form using Grassman-valued bond variables, and show insensitivity of the ground state 
density matrix to boundary conditions. 
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The distinction between different phases of matter is 
one of the most basic ideas in condensed matter and sta- 
tistical physics. In general, following Landau, to go be- 
tween two states of different symmetry in a classical sys- 
tem requires a phase transition, unless a symmetry break- 
ing field is turned on. In the case of quantum systems, 
similar ideas about symmetry hold. Suppose there ex- 
ist two local Hamiltonians, Tip and Tii, which both have 
an energy gap, AE, between a sector of approximately 
degenerate ground states and the rest of the spectrum. 
Then we ask: can one find a family of local Hamiltoni- 
ans, Hs, which depend continuously on s and interpolate 
between Tio and Ti.i such that the gap remains open? 

If so, then it is possible to go from Ho to Hi with- 
out a quantum phase transition. However, if so, then 
the technique of quasi- adiabaticQ continuation can be 
applied and leads, in certain cases, to a proof that there 
must be a phase transition due to an obstruction. For 
example in a ferromagnetic transverse field Ising model, 
^ = JT,^,Jn.n.SlSi + BJ2^Sl, One cau show that it 
is not possible to continue from a Hamiltonian with 
J >> B in the ferromagnetic phase to one with B >> J 
in the paramagnetic phase without a phase transition un- 
less one violates the Ising symmetry: one cannot have all 
of the Hs invariant under a global flip of all the spins. 
This is a case of a symmetry obstruction. 

However, there can be even worse obstructions. In the 
case of a topological obstruction, such as in the fractional 
quantum Hall effect, where Ho has a multiply degenerate 
ground state on a torus, any path in parameter Hs, if 
it includes only local Hamiltonians, can only produce an 
exponentially small splitting in the ground state sector [J, 
[2| unless the gap to the rest of the spectrum is closed. 
Here, symmetry does not play a role. A generalization 
of this kind of topological obstruction was presented in 
0, where it was shown that if all local operators have 
exponentially small matrix elements between the ground 
states for Ho, then this is preserved under quasi-adiabatic 
continuation. 

In general, we will take the ability to interpolate be- 
tween two different Hamiltonians, while preserving ap- 



propriate symmetries, without closing the gap as a def- 
inition of what it means to be in the same phase. Our 
main result here is to construct such a path in parameter 
space for a simple but important class of system: free 
fermion systems. These systems in some cases possess 
topological obstructions to continuation back to a triv- 
ial system, where we call a system trivial if there are no 
terms in the Hamiltonian coupling different sites. Exam- 
ples of such topological obstructions include Major ana 
number and Chern number The key idea of the con- 
struction in this paper is to evade these obstructions by 
doubling the system, taking two copies which are time- 
reverses of each other. 

Then, given the path in parameter space Hs, for s = 
to 1, between the original Hamiltonian Ho and a trivial 
Hamiltonian Hi, we use quasi-adiabatic continuation to 
write the ground state of the original system in a ma- 
trix product form[l| and show insensitivity to boundary 
conditions, a question naturally raised by work in and 
[7| where similar results were found for bosonic systems; 
another, less-efficient, way to write the ground state in a 
matrix product form is using thermal operators following 
[13] ■ Finally, in the conclusion, we discuss the relation 
between these results and the existence of exponentially 
localized Wannier functions. 

Family of Hamiltonians — We start by some definitions 
of the lattice and the Hamiltonian, then present the con- 
tinuous family of Hamiltonians, and finally give the state- 
ment and then the proof of the first result. Consider a 
lattice of V sites, labeled by .... Let 

H^Y. + E ^^^^1*] + ^-c- (1) 

be the free fermionic lattice Hamiltonian. The Hilbert 
space on each site i has two possible states, and 
are the creation and annihilation operators. 7i is a many- 
body operator, while Hij represent the matrix elements 
of a Hermitian V-hy-V matrix H and Ay are the matrix 
elements of an antisymmetric V-hy-V matrix A. Let 
dist(i, j) be a metric on the lattice. Assume the spectrum 
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of Ti. has a gap /S.E between the ground state and the first 
excited state. Let A be the matrix 



A: 



H A 

At -H' 



so that 



(2) 



(3) 



where (^,^t) jg the 2V dimensional vector 
(^fi, *2, This implies that 

the eigenvalue of A which is smallest in absolute value is 
equal to the gap AE in absolute value. 

We assume that Ti. is local in the following 
sense: for some constant /i > 0, and for some 
constant Si, then for all i, we have the bound 

2E,exp[Mdist(z,j)]yp7i7pT^~F < si < oo-,-,™^ 
implies the existence of a Lieb- Robinson bound|8l.l9lllCl|: 
there is some velocity v and some constant c which de- 
pend only on /i, si such that for any operators OxiOy 
with support on sets X,Y we have |j[Ox(0: Oy]|| < 
cx |X|||Ox||||Oy||exp[-/idist(X,r)](exp(z;Ai|i|)- 1). 

We introduce a system defined on two copies of the 
original lattice, with the same two-state Hilbert space 
on each site of each copy. We label sites on this system 
by (j, a), where i has V possible values and a —] , J,, and 
*1 ai ^i,a. are the creation and annihilation operators. We 
use the metric dist'((i, a), (j, h)) = dist(i, j) -I- (1 — 5a,b)- 

We set Tio to be the sum of Ti. on the first copy and 
Ti* on the second copy. We now define the continuous 
family of Hamiltonians with this Hilbert space: 



(4) 



The results we now show are that: (1): for all s, 
< s < 1, the Hamiltonian Tig has a gap equal to 
A£'; (2): for all s, < s < 1, the Hamiltonian 
Tis obeys the same Lieb-Robinson bound as TL, namely 
2(VT^Ejexp[Mist(z, j)]^|i/,,|2 + |A,,f+s2^i;) < 
si < oo; (3): the ground state of Tii is a product wave- 
function of the form 



^V2" 



^))|0), 



V2 



(5) 



where |0) is the state where all sites are empty. The 
proof of (3) is immediate. To prove (2), note that A_E 
is bounded above by si so that we may estimate the 
operator norm of sAi?(^'j i^*! | + h.c). 

We finally consider (1). Define the matrix Cs by 



C, = 



^/l^A 
sAi;i 



sAi;i 
-VT^A 



(6) 



Then, 

Us - \{^^,^\,^\,^^)^ ■ Cs ■ (vI/T,^'!,*},*^). (7) 

For each eigenvalue A of A, Cs has two eigenvalues equal 
to ± — s2)A2 + ,s2A£'2. Since the smallest eigenvalue 
of A was equal to AE in absolute value, the same holds 
for Cg and hence (1) follows. 

Finally, we comment on the reason for doubling the 
system, presenting two examples of topological obstruc- 
tions. The first example is a one-dimensional example 
based on the idea of Majorana number We consider a 
periodic chain of V sites, labeled j = 1...V. We define the 
operators C2j-i = '^j + and C2j — {^j — ^'])/«, thus 
giving Majorana operators Cfc defined for k = 1...2V. A 
quadratic Hamiltonian for a Majorana system is defined 
by a matrix A as in Eq. ([2]), while the ground state \l/o 
of the Hamiltonian has correlators given by 



(8) 



where the antisymmetric matrix B obeying = — 1 
is given by B = —isgn{A). Consider the state "if odd 
defined by the matrix Bodd with {Bodd)ij = S^+ij if 



1 mod 2 and {Bodd)i 



if J 



1 mod 2. 



where 1 is the unit 2y-by-2y matrix. 



This state is the state in which every site is occupied: 
(*odd|*I*j|^'odd) = 1- Now consider the state Beven de- 
fined by {Beven)ij = iBodd)i+i,j+i- Kitacv defines a Ma- 
jorana number for any anti-symmetric matrix B obeying 
= — 1 such that Bjk decays sufficiently rapidly in 
dist(j, fc). The Majorana number is an integer equal to 
±1 for an infinite chain, and is close to ±1 for a finite 
chain depending on how the chain length compares to 
the rate at which the coefficients decay. 

The Majorana number has opposite signs for the states 
^odd, ^even- This means that there is no way to find a 
family of free fermion Hamiltonians Hs which are local 
and gapped which have ground state ^even for s = and 
^odd for s = 1. To show this, suppose that such a family 
did exist. Then, given a gapped Hamiltonian, the matrix 
Bs has matrix elements {Bs)ij decaying exponentially 
rapidly in \ due to locality of correlation functions [If | 
and hence has a well-defined Majorana number up to cor- 
rections which tend to zero as V becomes large compared 
to the correlation length, and hence the Majorana num- 
ber cannot change sign in this evolution. 

However, doubling the chain evades this topological 
obstruction as the Majorana number is always even for 
the doubled system. The second example is in two 
dimension [5| • Consider a matrix P such that = P and 
such that P'^f^ is sufficiently rapidly decaying in dist(j, k). 
Then Kitaev defines a quantity v^P) for such a matrix 
which generalizes the notion of Chern number. This 
quantity i'{P) is shown to be equal to an integer in an 
infinite system, and to be close to an integer in a finite 
system, the error depending on how the decay of the co- 
efficients in dist(j, k) compares to the system size. 
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For any gapped local Hamiltonian of the form Ti, = 
'^\Hij'^j (so that there are no pairing terms: A = 0), 
we can define a projector Pjk — (^ol^'j^'fcl^'o)- Since 
the Hamiltonian is gapped, the correlations in the ground 
state ^0 will be short range, and hence Pjk will decay ex- 
ponentially rapidly in dist(j, k). Hence, if the system size 
is sufficiently large compared to the correlation length, 
v{P) will be close to an integer. A factorized wavefunc- 
tion has v{P) = 0, so therefore there is no way to con- 
tinue from a gapped Hamiltonians with non-vanishing 
v{P) to a Hamiltonian with a factorized ground state 
without either closing the gap or violating locality. How- 
ever, by doubling the system, and using the complex 
Hamiltonian for the second of the two copies, we evade 
this obstruction, since the total Chern number of the P 
defined by the ground state of Tis system vanishes for all 
s: at s = it cancels between the two copies. 

Quasi- Adiabatic Continuation and Matrix Product 
Ground State — The existence of a family of Hamilto- 
nians ([4]) implies that the ground state of Ho can be 
represented as a matrix product state as follows. The 
results that follow are valid for arbitrary families of lo- 
cal Hamiltonians Tig which have a gap for all < s < 1 
and which have a trivial ground state at s = 1, so we 
develop them in generality. Let / be some arbitrary 
length scale. Let us write an arbitrary local Hamilto- 
nian Tig as J^z'^zis), where Hzis) has support only 
on set Z; in this specific problem of free fermions, Z 
will always be a set of one or two sites. Using the tech- 
nique of quasi-adiabatic continuation, we can define a 
Hermitian operator T) such that ds'^o{s) — iVzis) + 
E^O(|^|exp[-^/e]^/Vll|5«^f^(s)||/A^;), where the 
correlation length ^ = 2v/ AE + /i. The operator I?z(s) 
has support on the set of sites within distance / of Z, 
and has \\Vz{s)\\ < ^/T/l(\\^sHzis)\\/ AE). For a free 
fermion system, as here, the operator is a bilinear 

in the fermion operators. 

Then, for sufficiently large Z, we can get a good ap- 
proximation to the 5'o(0) by 

^-0(0) « exp[z J Vs,ds']^o{l), (9) 

and then we can approximate the exponential by a ma- 
trix product operator of some given bond dimension, 
so that we can approximate ^'o as a matrix product 
state. It remains only to estimate the errors involved 
in the particular case to determine how large the bond 
dimension of the matrix product state must be. For 
(HI), we have \Z\ < 2 and \\dsHz{s)\\ < si, so for some 
ci which is a numeric constant of order unity, then for 
I = ci log(V^si/Ai?)^ + C2$, Eq. ([9]) gives an approxima- 
tion to ^'o(O) with error of order exp[— C2]. The unitary 
operator exp[z J^Vs'ds'] can be regarded as time evolu- 
tion under an effective time dependent Hamiltonian Vgr . 
This effective Hamiltonian has a Lieb- Robinson group ve- 
locity which we denote ^' which is bounded by a constant 



of order unity times sil/AE. Then, for the given I, we 
can approximate exp[i I?s'ds']5'i(s) within error e by 
a matrix product operator Umps as follows. 

This kind of approximation was discussed for one di- 
mensional spin systems in [l3| . The construction of [l3l | 
proceeds by breaking the one dimensional system into 
short intervals labeled 1,2,..., and showing that the uni- 
tary operator can be approximately written as a quantum 
circuit: Vi,2V2.3-..UiU2--.., where the operators Ui, U2, ... 
act only on each interval 1,2, ... and the unitary operators 
Vi^i+i act only on the right-half of interval i and the left 
half of interval i + 1, so that [Ui,Uj] = [Vi^i+i, V^.j+i] = 0. 
This construction can be directly generalized to higher 
dimensions [l3| • As a next step, auxiliary bond variables 
are introduced to write the operators V^^i+i as a sum of 
product of operators on interval i and interval i + 

In our problem, in order to decompose the operator 
Vi^i+i into a sum of products of operators, we must in- 
troduce Grassman-valued bond variables. In order to 
approximation the unitary evolution to an error e, we 
need the number of Grassman-valued bond variables 
on each bond to be logarithmically large in V, 1 /e and 
proportional to ^' (for the spin system, the bond di- 
mension is exponentially large in this quantity). Thus, 
we can approximate ^o(O) by a matrix product state 
^mps — C/mps ^o(l) with Grassman-valued bond vari- 
ables. We note that the bond dimension required for 
this free fermion system is much smaller than that for 
quasi-adiabatic continuation of a general interacting sys- 
tem, where a bond dimension polylogarithmically large 
in V may be required, because "Dz is a fermion bilinear. 

Sensitivity to Boundary Conditions — The final ques- 
tion we consider is the sensitivity of such a gapped 
fermion system (j4|) to boundary conditions at s = 0. 
Consider two systems with Hamiltonians TL,TL' which 
agree on some set of sites X in the following sense: for any 
operator O with support on X, [0,H] = [0,W]. Now, 
let Y be some subset of X such that the set of all points i 
with dist(z, y) < Hs a subset of AT. We then consider the 
difference between the reduced density matrices, py,p'yi 
of the two systems. We will show that as I gets large, the 
difference between the two density matrices converges to 
zero in trace norm. Physically, this may be viewed as 
a statement about insensitivity to boundary conditions: 
even if two systems differ at the "boundary" (that is, 
outside X), far enough away from the boundary (that is, 
inside Y for large enough I) the physical properties will 
be the same. 

We again use quasi-adiabatic continuation. We take 
two copies of each system, and define a continuous fam- 
ily of Hamiltonians HsjTi's as above, so that Hi = H[. 
Let 4'o(s),^n(s) be the ground states of these Hamilto- 
nians. Then[]| for any operator O with support on Y we 
can define an operator 0{s) such that 0{s) is supported 
on X |(*o(l)|0(l)|*o(l)) - («'o(0)|0|*o(0))| < e where 
the error e is exponentially small in Since Hs,H'g 
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agree on X, we find that also \{%{1)\0{1)\%{1)) - 
{%{0)\O\%{0))\ < e. However, since Hi = n[, we 
have («'o(l)|0(l)|4'o(l)) = (^(1)10(1)1*^,(1)), and so 
|(*o(0)|0|*o(0)) - {%mO\%{0))\ < 2e for any O with 
support on Y. Therefore, we can bound the trace norm 
distance between pY,i and py,2- ^^{\py,i — Py,2\) < 2e. 
We can estimate e in the specific case of a d-dimensional 
free fermionic system to get 



tr(|py,i~py,2|) < const. X \Y\^l/vAEsi exp[-//e'](^')''- 

(10) 

Discussion and Wannier Functions — Can there be 
topological obstructions to continuing a time-reversal 
symmetric undoubled system back to a trivial system? 
One possible example is the spin Hall state [l^, but by 
turning on an external magnetic field, breaking the time 
reversal symmetry, it is possible to continuously deform 
the undoubled Hamiltonian to a trivial Hamiltonian. 

In order to understand possible obstructions in time 
reversal symmetric systems better, we consider Wan- 
nier functions. We say that a system has exponentially 
localized Wannier functions if one can find a basis of 
states \(j)n) which span the occupied states of the free 
fermion system (so that the Green's function is equal to 
G = ^„ \ijjn) {(t>n\) and which arc exponentially localized 
in real space. Such functions have long been considered 
for periodic systems, as well as more recently for disor- 
dered systems [17]. 

If a free fermion Hamiltonian Tio can be continued to 
a trivial Hamiltonian Tii, then exponentially localized 
Wannier functions exist: the Hamiltonian T-Ci has Wan- 
nier functions \4>n), where each |0„) has support on a sin- 
gle site. Then, by quasi-adiabatically continuing from Tii 
to Tio, we can find exponentially localized Wannier func- 
tions for Tio- Only very recently [l6| was it shown that 
for time reversal symmetric, periodic systems in two and 
three dimensions exponentially localized Wannier func- 
tions do exist, in that case without any doubling required. 
On the other hand, the proof in [16] suggests that there 
exist time-reversal symmetric, four-dimensional, periodic 
systems for which exponentially localized Wannier func- 
tions do not exist. These systems would be in a different 
topological phase from a trivial Hamiltonian, and even 
by adding a magnetic field it is not possible to continue 
back to the trivial Hamiltonian without closing a gap. 

An interesting open question is when such exponen- 
tially localized Wannier functions must exist for gapped, 
undoubled, time-reversal symmetric, disordered systems. 
In one dimension, with open boundary conditions, this 
can be shown as follows: let X be the position operator. 
Then, GXG is a Hermitian operator and can be diag- 
onalized; because G is short-range due to the gap [111], 
the eigenvectors of GXG are exponentially localized in 
real-space and hence can be used as Wannier functions. 
In two dimensions, we consider the pair of operators 



GXG, GYG. These operators almost commute, since 
\\[GXG,GYG]\\ ~ 0(1) « \\GXG\\ ~ L, where L is 
the linear size of the system. If these operators are are 
close to a pair of commuting operators we could simulta- 
neously diagonalize those operators and use their eigen- 
functions as Wannier functions. The question of when 
almost commuting operators are close to commuting op- 
erators is well-studied [13], and those results can be used 
to show that given a family of systems of increasing size 
L and constant gap, then for any 5 > for sufficiently 
large L we can define Wannier functions localized on a 
scale SL. In three dimension, we have three almost com- 
muting operators GXG, GYG, GZG, which need not be 
close to commuting operators as shown in JT^. 
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